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We study the development of mean structures in a non¬ 
linear model of large scale ocean dynamics with bottom to¬ 
pography and dissipation, and forced with a noise term. We 
show that the presence of noise in this nonlinear model leads 
to persistent average currents directed along isobaths. At 
variance with previous works we use a scale unselective dissi¬ 
pation, so that the phenomenon can not be explained in terms 
of minimum enstrophy states. The effect requires the pres¬ 
ence of both the nonlinear and the random terms, and can be 
though of as an ordering of the stochastic energy input by the 
combined effect of nonlinearity and topography. The statisti¬ 
cally steady state is well described by a generalized canonical 
equilibrium with mean energy and enstrophy determined by a 
balance between random forcing and dissipation. This result 
allows predicting the strengh of the noise-sustained currents. 
Finally we discuss the relevance that these noise-induced cur¬ 
rents could have on real ocean circulation. 

PACS numbers: 05.40.+j, 92.90.+X, 47.27.Eq 


I. INTRODUCTION 


Under certain conditions, nonlinear interactions are 
able to organize random inputs of energy into directed 
motions |i],[|. These ordered fluxes or structures sup¬ 
ported by noise have important effects on specific prob¬ 
lems in biology ||[|], engineering [g or physics @.[7J. In 
fluid dynamics noise appears in a variety of contexts. 
Thermal fluctuations give rise to Langevin source terms 
in the equations of motion g|J , while stirring forces can 
be modeled as stochastic terms forcing the fluid motion 
f[0[ . Noise terms have also been used in modeling short- 
scale instabilities such as the baroclinic instability in geo¬ 
physical problems 00 - A different way in which noise 
can enter the description of fluid systems is by consid¬ 
ering the interactions taking place at scales below the 
resolution of a computer model. Addition of noise is a 
(crude) way of including some influence from unresolved 
scales [ 00 - 

The influence of small noise on large scale fluid dy¬ 
namics, far enough from hydrodynamic instabilities and 
phase transitions, is expected to be small. However we 
address in this paper a situation in which it has impor¬ 
tant consequences: the coupling of fluid to bottom topog¬ 
raphy in ocean dynamics, which is one of the most im¬ 
portant but poorly understood factors controlling ocean 


circulation. We will show that, when interacting with 
noise, this coupling leads to the development of steady 
currents, even in the absence of steady forcing. Although 
our results are derived for a particular quasigeostrophic 
fluid model, we expect our general conclusions to apply to 
other problems in which fluid motion occurs on irregular 
geometries in the presence of noise. 

As a first step the work in this paper will be restricted 
to barotropic, randomly forced, two-dimensional quasi¬ 
geostrophic turbulence. Several reasons have been con¬ 
sidered for this selection. First of all, this framework pro¬ 
vides a physical model complex enough to contain some 
relevant aspects of ocean circulation, but easy enough for 
allowing physical understanding without the need of too 
complex numerical experiments (always difficult to in¬ 
terpret). Second, baroclinic nature of ocean circulation 
can be usefully represented for large scales by randomly 
forced barotropic models |ll|. In these barotropic mod¬ 
els, the random forcing formulation is a simple Markovian 
process that provides a reasonable representation of the 
baroclinic nonlinear interactions fi~^JT^ . 

A brief summary of known results relevant to our study 
follows for completeness. In the last decades much ef¬ 
fort has been devoted towards understanding the effects 
of bottom topography in barotropic models. Most rel¬ 
evant results on this subject can be summarized in the 
several key works starting by that developed by Salmon 
et al. 0' where the statistical properties of topographic 
turbulence neglecting dissipation and forcing were stud¬ 
ied. Among their results the demonstration of the ten¬ 
dency to reach a maximum entropy (Gibbs) state char¬ 
acterized by the existence of stationary mean currents 
following isobaths should be stressed. The role of dissi¬ 
pation was analyzed in Jl6[ using a different approach. In 
their study viscous damping was hypothesized, due to its 
selective damping of high wavenumber modes, to drive 
the flow towards a state of minimum potential enstro¬ 
phy for a given energy, hence dissipating enstrophy more 
rapidly than energy which remains almost constant for a 
long time. This long-lived minimum potential enstrophy 
state also consists of flow along topographic contours. 
More recently, stability properties of minimum enstro¬ 
phy and maximum entropy states and their relationship 
have been analyzed by Carnevale and Frederiksen [|17|. 
A quantitative theory of time-dependent, forced and dis¬ 
sipative flow over topography was analyzed by Herring 
[p~8| who followed the direct interaction approximation 
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(DIA) [Q and a modified form of the test-field-model 
turbulent formalisms J|(J. He found that in the case of 
random stirring and viscous topographic turbulence the 
stationary components of the energy spectrum behaved 
close to the above mentioned inviscicl results of Salmon 
et al. [ |l5|. A statistical theory simpler than DIA was 
used in |2l[| to investigate the time evolution of the en¬ 
ergy and vorticity-topography correlation spectrum. His 
results also show the emergence of steady circulation for 
a certain range of the parameters of the theory. As a 
first step towards understanding the more complicated 
baroclinic case, Treguier j|2^| made an extensive explo¬ 
ration of the parameter space for barotropic turbulence. 
Her numerical results show a non-monotonic dependence 
of the percentage of steady currents on the forcing, vis¬ 
cous and topography parameters. In some simulations 
where random forcing and dissipation were considered, 
these steady currents resemble the simple inviscid solu¬ 
tion found in fL5| . 

All these previous studies address the existence of sta¬ 
tionary currents following isobaths in a variety of condi¬ 
tions involving random noise. However from these studies 
it is not possible to infer a clear relation between noise 
and the physical nature of these currents. The reason is 
that scale selective (i.e., viscosity) dampings were con¬ 
sidered so that the generation of directed currents can 
be explained in terms of minimum enstrophy states [ fTof . 
This physical mechanism involving minimum enstrophy 
solutions does not require the existence of noise. 

In this paper, we try to highlight the role of noise in 
the appearance of steady currents in such open systems. 
The main difference with previous studies is that we have 
considered topographic turbulence damped by Rayleigh 
friction (bottom friction) instead of viscosity damping. 
This introduces important physical differences because 
Rayleigh friction provides a scale-unselective way of dis¬ 
sipating energy from the system and thus does not pro¬ 
duce transient minimum enstrophy states. In this way, 
avoiding one of the possible mechanisms for explaining 
the mean flows that were present in the simulations of 
[PI , and [p2|, we show that the combined effect of 
noise and nonlinear terms is enough to generate mean 
flows when bottom topography is present. In order to 
get a more physical understanding of this phenomenon, 
we show that this current formation at large scales is de¬ 
scribed, to a good degree of approximation, by a general¬ 
ized canonical distribution in which the mean enstrophy 
and the mean energy are fixed by the balance between 
random forcing and dissipation. Then we suggest that 
the generation of currents by noise in this system can 
be understood as arising from the tendency of the large 
scale fluid motions in our model to approach equilibrium¬ 
like states, that is states with the maximum possible en¬ 
tropy consistent with given constraints. However, the 
non-equilibrium character of our system is reflected in 
the small scales, which depart from this equilibrium-like 
state. 

Two comments are in order: First, since our motiva¬ 


tion is the understanding of large scale ocean circulation, 
we are interested in the final asymptotic behavior of the 
fluid system under forcing and dissipation, and not in the 
approach to it. Therefore the numerous studies available 
for freely decaying turbulence can not be applied to this 
problem. Second, the presence of topography introduces 
spatial scales into the problem so that mechanisms re¬ 
lated to scale-invariant wavenumber cascades [|23f can not 
be as relevant here as in a flat bottom situation. 

The paper is organized as follows: in Section |h| the 
quasigeostrophic model is introduced and numerically 
solved to show the appearance of directed currents along 
the topography. Section III contains analytical consid¬ 
erations showing the relevance of the nonlinear terms, 
discussion of the Fokker-Planck equation for the model 
and calculation of some average values. The generalized 
canonical distribution is also presented and some conse¬ 
quences are derived from it. Section IV checks the valid¬ 
ity of the canonical description. Discussion of the results 
and their implications for large-scale ocean modelling are 
in Section 
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II. THE MODEL AND THE NOISE-SUSTAINED 
CURRENTS 


The quasi-geostrophic evolution equation for the 
streamfunction i/j(x, y , t) describing flow over topography 
is given by [ pdf : 

^ + ,7(^VV + h') =-eV 2 <p + F, (1) 

where 

J (A, B) = k • (VA x VB) (2) 

e” 1 is the bottom friction decay time, h! = fAH/Ho 
with / the Coriolis parameter (that will take the value 
10~ 4 s _1 as appropriate for mean latitudes), Hq the mean 
depth, and AH(x,y) the local deviation from this mean 
depth, k is the unit vector normal to the horizontal plane 
{x,y), x denotes the cross vector product and F(x,y,t) 
is a Gaussian white-noise process with zero mean and 
correlations (F(x, y, t)F(x', y\ t')) = DS(x — x')S(y — 
y')5{t — t'). It represents relative-vorticity random forc¬ 
ing. Eq.(m) is to be solved in a square domain with 
boundary conditions such that if} vanishes at the bound¬ 
aries. The streamfunction il> gives the horizontal compo¬ 
nents of the fluid velocity (u(x,y),v(x,y)) as 


dij) d'i/j 

dy dx 


( 3 ) 


In order to solve (jlj) we have used the quasi-geostrophic 
numerical scheme developed in |2^| on a grid of 128 x 128 
points. The resolution of our scheme (corresponding to 
the distance between grid points) is of 10 km, so that 
the total system size is L = 1280 km. The algorithm, 
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based on finite differences and the leap-frog algorithm, 
keeps the value of energy and enstrophy constant when 
it is run in the inviscid and unforced case. The consistent 
way of introducing the stochastic term into the leap-frog 
scheme can be found following the lines of |2^,|2^] : Basi¬ 
cally, a field of Gaussian random numbers is generated at 
each time step, and the random field added to the deter¬ 
ministically evolved vorticity field at each time step is the 
sum of the field generated at this time step plus the field 
generated at the previous time step. The random num¬ 
ber generator employed in the simulation is detailed in 
p8f . The amplitude of the forcing, D = 2 x 10 -9 m 2 /s 2 , 
and the e-folding decay constant, e = 3 x 10 -7 s -1 of the 
damping mechanism have been chosen in order to obtain 
final velocities of several centimeters per second. The 
topographic field (Fig 1) is randomly generated from a 
isotropic spectrum containing, with equal amplitude and 
random phases, all the Fourier modes corresponding to 
scales between 80 km and 300 km. The model was run 
for 5 x 10 5 time steps (corresponding to 126 years) after 
a statistically stationary state was reached. The average 
streamfunction was then computed from the configura¬ 
tions obtained every time step. 

The results obtained from this numerical experiment 
show the generation of a stationary mean flow correlated 
with bottom topography and with velocity values of or¬ 
der of 3 cm/s (Fig. 2). This numerical result indicates 
the existence of a physical mechanism, different from the 
minimum enstrophy tendency addressed in , that is 
able to organize random inputs of energy in a mean state 
different from rest characterized by currents with a clear 
dependence with bottom topography. The mean currents 
will advect any neutral tracer with them, so that the kind 
of motion so generated is very similar to the phenom¬ 
ena described in 01- In addition we have checked by 
Fourier analysis that the temporal series of ip(x, y, t ) at 
fixed space point ( x , y ) is roughly periodic in time (but 
quite noisy), so that there are also directed currents in 
the phase space of ip. In the next section we try to clarify 
the possible origin of the currents. 


k • (V (ip) x Vti) + eV 2 (ip) = 0 


( 5 ) 


No stationary solution of Eq. ((H) can be a function of 
topography. This is because a functional dependence of 
the mean streamfunction on the bottom topography im¬ 
plies that the first term of (|B|) would be identically zero 
and then < ip > should be a solution of Laplace’s equa¬ 
tion. With the boundary conditions used, this implies 
that the mean streamfunction vanishes in every point of 
the system. Thus in the linear model (^) the forcing gen¬ 
erates small excursions of < ip > around the rest state, in 
a way unrelated to topography. In conclusion, the nature 
of noise sustained currents is thus determined not only by 
the random noise but also by the nonlinear interactions. 


B. The Fokker-Planck equation of the nonlinear 
model 

As it was shown from the study of the linearized dy¬ 
namics, the nonlinear term is needed in order to get mean 
stationary solutions different from rest. This fact im¬ 
plies to consider the full dynamics given by ([!]). Due to 
the stochastic nature of the system we will explore the 
statistics of its solutions in phase space. To do that, the 
Fokker-Planck equation (FPE) of the probability distri¬ 
bution of the field solution of Eq. (jlj) will be derived. 
From the standpoint of treating the statistics of an en¬ 
semble of solutions of (Q), it is convenient to represent 
the streamfunction ip(x,y,t) as a linear combination of 
orthogonal eigenfunctions of the Laplacian satisfying the 
adequate boundary conditions: 

ip(x,y,t) = ^Ai(t)(pi{x,y), V 2 cpi(x,y) = - a 2 (pi(x,y) 

i 

( 6 ) 


dxdy<pi(x, y)<pj (x, y) = S6 i: 


( 7 ) 


III. ANALYTIC APPROACHES 
A. The linear model 

As a first step to investigate the emergence of these 
noise induced currents, we will explore the solutions of 
the linearized dynamics. Neglecting the nonlinear term, 
Eq. (|Ij) transforms into: 


d\7 2 ip 

dt 


+ k • (yip x V/i') = —tS7 2 ip + F 


( 4 ) 


After ensemble averaging Eq. (|j) over noise realizations 
we obtain the equation that determines the stationary 
state of the linear system: 


The integral is over the whole two-dimensional domain, 
of area S = L 2 . 

The evolution equation for the amplitude factors Aj (t) 
reads: 

a k A k -5> fc (aijAjAi - Aihj) = - ea 2 k A k + a k f k 


where: 




D 


Sa k a k ‘ 


■S k k'8{t - t') , 


( 8 ) 


akfk = —^ J dxdy F(x, y,t)<Pk(x,y) , (9) 


( 10 ) 
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and 


Pijk = J dxdy k- (V(/>j(x, y) x V(j)j(x, y)) <f> k (x, y) . 

( 11 ) 

Some properties of the nonlinear interaction coeffi¬ 
cients flijk are that they vanish if any two indexes are 
equal, they are invariant under cyclic permutation of in¬ 
dices, and they reverse sign under non-cyclic permutation 
of indices |f29| . 

In order to use the symmetries and asymmetries of 
Eq. (||) more clearly, we introduce the following change 
of variables: 


X k (t) = a k A k (t) , hj = cxjhj . ( 12 ) 

The equation is now: 

X k (t) = N k - eX k (t) + f k (t) (13) 

where: 

Nk=Y (.XiXj - Xihj) (14) 

OtjCtiOtk 

IJ 

Now, we consider the amplitudes {W} of the orthogo¬ 
nal modes as the coordinates of a point (in an infinite¬ 
dimensional phase space) whose position corresponds to 
the state of a particular realization at a single instant 
in time. Since f k (t) is white noise (Eq. ©), the con¬ 
tinuity equation for the phase space probability density 
p = p\{Xi},t) takes the Fokker-Planck form [30| : 

m + ^~dT k p ( Nk ~ eXk) = (15) 

k k K 


where p k = D/(2Sa k ) and represents a diffusion coeffi¬ 
cient for the A:—mode. We can see from @) that random 
forcing acts as a diffusion mechanism for the probability 
density in phase space. This mixing character is oppo¬ 
site to the term represented by friction that tends to drive 
the system towards the origin of the phase space, the rest 
state. 

The FPE (|l5|) can be easily integrated in the stationary 
state if the potential condition: 


dZ k dZi 
dX, 8X k 


(16) 


with Z k = p^ 1 (N k — eX k ) holds |j(J. Otherwise there is 
no simple way of finding the stationary probability dis¬ 
tribution. The potential condition is not satisfied in the 
case described by m so that some assumptions must be 
performed to get an approximate probability distribution 
for the stationary state. 

However the FPE can be still useful for determining 
some average quantities needed below. The average of 
any phase space function F({X k }) satisfy the adjoint of 


Eq.(|l5|), or backwards equation Q. This can be shown 
by multiplying (|l5| ) by F({Xfc}), integrating over phase 
space, and using integration by parts. For stationary 
averages (|j = 0), the result is 

r / f)F fP F \ 

J dT^—( N k -eX k )+p k -^)p = 0 (17) 

f dr denotes integration over phase space JYl k dX k . 

Choosing F = \a k (\X k — h . k )~, substituting into (|l7| ) 
and summing over k. the antisymmetric terms containing 
f3ij k disappear and we find: 

e J2 a k((\ Xk ~h k ] X k \ =J2^p k al (18) 

where < ... > indicates average over the stationary dis¬ 
tribution p. The term between brackets is related to the 
mean potential enstrophy, so that Eq. ( |T8| ) gives an ex¬ 
pression for this quantity in terms of the damping coeffi¬ 
cient e and the noise intensity p k ■ Similarly, by choosing 
F = \X k and repeating the procedure above we find: 

e J2 a ‘ k (\ x k') = H\dk (i9) 

k ' k 

This relates the mean energy {\X^) to e and p k . 


C. A generalized canonical equilibrium 

Properties of the stationary state of the FPE (jl5|) can 
be used to obtain an approximate expression for the prob¬ 
ability distribution. As it was shown in jn|, invariance 
of energy and potential enstrophy has important impli¬ 
cations on the inviscid and unforced dynamics. In par¬ 
ticular statistical properties can be correctly computed 
by assuming that at long times the system explores com¬ 
pletely, or at least samples with enough significance, the 
hypersurface of constant energy and potential enstro¬ 
phy (for a discussion see |l2)). When damping is intro¬ 
duced, the system will be brought to the vicinity of the 
rest state in phase space. But random forcing will keep 
some motion present. Since the effect of the nonlinear 
terms in the unforced and undamped dynamics was to 
drive the system along the intersection of the constant- 
energy and constant-potential-enstrophy hypersurfaces, 
it is then reasonable to expect that in the presence of 
forcing and damping evolution will proceed on average 
close to the hypersurfaces of mean constant energy and 
mean potential enstrophy determined by the damping 
and forcing coefficients through equations (|l^ ) and (Hj). 
If this hypothesis holds, the most likely probability dis¬ 
tribution compatible with such constraints, that is the 
one obtained from the maximization of the entropy func¬ 
tional S[p] = J Y[ k dX k p{\X k }) logp({X fc }) subjected to 
the constraints ([l8]) and should be an approxima¬ 
tion to the real stationary probability distribution 
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(25) 


The validity of the heuristic considerations above, as 
the validity of general statistical mechanics arguments 
in any interesting physical situation, would be hard to 
establish rigorously. In our case, justification of our as¬ 
sumptions will come from comparisons with numerical 
simulations. Ergodicity, or at least efficient sampling of 
the constant energy and potential enstrophy hypersur¬ 
faces would be certainly more justified here than in the 
unforced case because of the additional mixing tendency 
provided by random forcing. 

Entropy maximization subjected to constraints (|I|) 
and © leads formally to the same expression as in the 
inviscid and unforced case, that is, the predicted proba¬ 
bility distribution is: 

p({X k }) = l[p k (X k ), (20) 

k 


with 

p k (X k ) = e {-h{ a + ba l) x l+ a l h kb x k) (21) 


a and b are the Lagrange multipliers enforcing the con¬ 
straints. Although this expression is formally identical 
to the inviscid one, there is an important difference: In 
the inviscid case the Lagrange multipliers are determined 
from the initial values of energy and enstrophy, whereas 
here the dependence on the initial state is asymptoti¬ 
cally lost and (||l|) refers to the final state in which the 
Lagrange multipliers are determined by the mean values 
of energy and enstrophy. Their exact values are given in 
(|l|) and ( p~9| ) in terms of the forcing and friction param¬ 
eters present in the original equation (|l|). The values of 
the first and second moments of X k obtained from ( |2l| ) 
are given by: 


(*k) = 


alh k b 

a + ba 2 k 


( 22 ) 


{Xl) = 


a 


ba l 


ba l 


(23) 


These equations in conjunction with (|18| ) and (|19|) al¬ 
low the determination of the Lagrange multipliers a and b 
from the parameters in the FPE or in (ll|). Alternatively, 
they relate the Lagrange multipliers with quantities mea¬ 
surable from simulation or experiments. 


As in the inviscid case, equations (22) and 


ply the existence of a mean circulation controlled by 
the topography |l5|]. To see that we focus on the value 
of the first moment which is specially interesting since 
it describes the generation of mean flows following iso¬ 
baths. This becomes evident rewriting ( p2|) in terms of 
the Fourier amplitudes of the functional expansion (H): 


( A k) = 


h'b 


ba l 


(24) 


that in the physical space imply the relation: 


V 2 (ip) +h' = c (ip) 

where c = a/b. This linear relationship between the mean 
potential vorticity and the mean streamfunction implies, 
considering that for large scale motions the relative vor¬ 
ticity is smaller than the ambient one (V 2 (ip) <C h '), that 
(ip) is proportional to h!, that is (ip) describes mean flows 
followingthe isobaths. We have checked that the distri¬ 
bution (go|)-(^lj) is not an exact solution of the FPE dH)- 
In the following Section however we show that it gives a 
good description of the large scale average properties. 


IV. COMPARISON WITH NUMERICAL 
SIMULATIONS 

The correlation between topography and mean stream- 
function, as predicted by (||^) is obvious by comparing 
Figs. 1 and 2. The agreement can be made more quanti¬ 
tative by plotting V 2 (ip) + h' versus (ip) , both obtained 
from the numerical simulation. This is done in Fig 3. 
Eq. (^5|) predicts a linear relationship. Fig. 3 shows 
a clear linear trend but there is also an evident scat¬ 
ter in the data that will be commented below. The lin¬ 
ear correlation coefficient is 0.940. Further analysis of 
the validity of theoretical predictions can be done by fit¬ 
ting the data in Fig. 3 to Eq. (||[) to get a numerical 
value for the constant c. This procedure gives a value of 
c = (5.30±0.0T) x 10~ 9 m~ 2 . The predicted value of c ob¬ 
tained from ( |22| ) and (^3|) by considering the magnitudes 
of the mean energy and potential enstrophy obtained in 
the stationary state is 5.5 x 10 _9 to^ 2 which represents 
a deviation of 4% with respect to the value obtained by 
numerical simulation. In order to test if such deviations 
occur at all scales, we have computed the power spec¬ 
trum of the fields (V 2 (ip) + h')/c and (ip), both obtained 
from the numerical model, and compared them with the 
energy spectrum of the field (ip) computed from direct 
resolution of ( p5| ) which will be taken as reference spec¬ 
trum. Theory predicts that these distributions should be 
the same. Fig. 4 shows that there is very good agreement 
between numerical results and the theoretical spectrum 
obtained from the canonical probability distribution at 
low wavenumbers (large scales) where topography com¬ 
ponents are relevant. The most important deviations oc¬ 
cur at small scales. There are also much smaller discrep¬ 
ancies at very low wavenumbers, that is at basin scales, 
similar to the ones observed when analyzing the equilib¬ 
rium solutions of purely inviscid two-dimensional turbu¬ 
lence J33| . The reason for them could be the long dynam¬ 
ical time scale of large-scale eddies that implies very long 
equilibration times [p3| . The deviations involving small 
scales cannot be explained by relaxation-time arguments. 
These scales are clearly out of equilibrium. This absence 
of equilibrium at least at small scales could be expected 
from arguments based on turbulent cascades: In two- 
dimensional turbulence an enstrophy flux towards small 
scales is expected Q. Since Rayleigh friction is not as 
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effective in dissipating vorticity at small scales as viscos¬ 
ity would be, an accumulation of enstrophy is expected at 
the largest wavenumbers, and this is in fact seen in Fig. 
4. In contrast, equations ©-© predict that states 
with vorticity at wavenumbers above a cut-off are for¬ 
bidden. This is easy to see from the canonical partition 
function of the system, that can be written as: 


Z{b a) = J Dip dxd y( v ^ 2 +i 2 f dxd v (y 2 ii>+h') 2 ) 

(26) 

where l 2 = b/a = c -1 has units of length square and 
f Dip denotes functional integration. If X/ 2 ip + h' varies 
on length scales smaller than l the second term of the 
exponential dominates and all the microstates of the sys¬ 
tem with these characteristics are dropped out from the 
ensemble averages. As a result properties of scales near 
or smaller than the cut-off cannot be well described by 
the canonical probability distribution. 

The fact that small scales are responsible for the broad¬ 
ening of the linear tendency in Fig. 3 can be verified 
from the fact that the linear correlation coefficient be¬ 
tween mean potential vorticity and mean streamfunction 
fields after scales smaller than 100 km (i.e. several times 
the cut-off scale) have been filtered out is 0.975, higher 
than the original one. Finally, if all scales except those 
with relevant topographic components (100 — 400 km) 
are removed, correlation coefficient increases until a value 
of 0.982, that is, the agreement with theory predictions 
is even better. In conclusion the canonical distribution 
given by (20)- © is a good description of the numerical 
results for length scales larger than 100 km. 


V. DISCUSSION 

In this paper, we have focussed our attention on the ef¬ 
fects that noise produces on a simplified stochastic quasi- 
geostrophic model where only friction and random terms 
have been considered. The reason for our interest comes 
from indications of generation of currents along isobaths 
in randomly forced quasigeostrophic systems presented 
in previous works |l^,^l|,^2| . Since the main interest of 
those works was not the study of the effect of random 
forcing on flow with underlying topography, the formu¬ 
lations used there included several effects absent in our 
model. In particular numerical experiments were carried 
out in the presence of correlated noise and scale-selective 
dampings. This more complex framework makes it diffi¬ 
cult to identify unambiguously the basic physical mech¬ 
anisms producing the currents. In particular the intro¬ 
duction of viscosity or bi-laplacian damping generates de¬ 
caying minimum enstrophy states described by ||l6| that 
obscure the effects of noise highlighted in our work. Our 
use of bottom friction as a damping mechanism tries to 
isolate the effect that noise could produce on the dynam¬ 
ics by avoiding minimum enstrophy states. 


As a result of our study, we have found that the ran¬ 
dom term supplies energy to the system and contributes 
(together with the nonlinear term) to the mixing of the 
probability density in phase space. This compensates the 
reducing volume tendency of the dissipation term. When 
this compensation occurs the system dynamics resembles 
the behaviour developed by the isolated system, that is 
without forcing and damping. The physical effect of this 
phase space dynamics is the generation of a stationary 
state characterized by mean fluid motions following iso¬ 
baths. We can conclude then that these mean currents 
are noise-sustained and come from the interaction of ran¬ 
dom noise with the nonlinear term. The existence of 
bottom topography is relevant in order to get a mean 
component in the spectrum, that is, to develop mean 
currents in the basin. This result seems to indicate that 
noise could play a relevant role in maintaining the mean 
component of planetary motions when bottom topogra¬ 
phy is present. A remaining open question is however, 
how much of the currents observed in [ ^8pl|j2^ | is pro¬ 
duced by random terms and how much is related to the 
minimum enstrophy tendency that is certainly present in 
their models. Our isolation of the effects of noise is a first 
step to evaluate the relative importance of both effects. 

One of the state of the art questions on the study of 
large scale ocean circulation and modeling is related to 
the effects that unresolved scales might have on the larger 
scale dynamics. The range of scales which numerical 
models can handle explicitly is limited by the actual com¬ 
puter power, the large size of the global domain and the 
time characteristics of climatic studies (long-time simu¬ 
lations). These limitations prevent us from simulating 
the small and large scale features simultaneously. How¬ 
ever, it is also well known that some of these unresolved 
processes produce significant influences on the large scale 
circulation. The influence of small scales must then be 
introduced through some parametrization. 

A natural way of considering the unresolved dynamics 
would be to include the small and fast oceanic processes 
as some kind of random noise in the deterministic mod¬ 
els of large scale phenomena. In general, dealing with 
ocean modeling the action of small scales is only consid¬ 
ered as an enhancement of the viscosity coefficient (the so 
called eddy-viscosity). Random noise is usually neglected 
from oceanic calculations because it is assumed that its 
presence produces a negligible effect on large scale ocean 
dynamics. 

Although the full ocean dynamics is characterized 
by the existence of three-dimensional motions, for suf¬ 
ficiently large spatial scales (scales larger than several 
hundreds of kilometers) and long time behaviour (more 
than one week) the ocean can be well described by quasi¬ 
geostrophic dynamics including bottom topography. The 
results shown here could be a first basis from which to 
conjecture that the origin of some observed large scale 
ocean structures are due to phenomena which are several 
orders of magnitude smaller than the structures them¬ 
selves. The experimental determination of this conjec- 
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ture can be very difficult to carry out due to the broad 
range of spatial and temporal scales in ocean. How¬ 
ever, new approaches for future studies where large scale 
stochastic primitive equation models are used, instead 
of purely deterministic ones, could be a reasonable way 
to go deeper into the problem. Effectively, inclusion of 
noise in a coarse-resolution primitive equation model as 
a representation of the small scale processes that are not 
well resolved would provide a mechanism for the model 
to generate noise-sustained large scale currents along iso¬ 
baths in a natural way, that obviously would be missing 
in a deterministic version of such a model. Ad-hoc inclu¬ 
sion of currents along isobaths has been already shown 
to improve results of coarse-resolution numerical mod¬ 
els Presently, general ocean circulation models 

are deterministic. However, based on the results estab¬ 
lished on this paper it seems reasonable to hypothesize 
that stochastic models will provide more realistic answers 
about ocean dynamics and more concretely about the 
possible existence and relevance of these noise-sustained 
currents in the real ocean. 
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FIG. 1. Depth contours of a randomly generated bottom 
topography. Maximum depth is 381.8m and minimum depth 
-381.8m over an average depth of 5000m. Levels are plotted 
every 63.6 m. Continuous contours are for positive deviations 
with respect to the mean, whereas dashed contours are for 
negative ones. 

FIG. 2. Mean streamfunction computed by time averag¬ 
ing when a stationary state has been achieved. Maximum 
and minimum values are 1120.2 and -1120.2 m 2 /s. Levels are 
plotted every 186.7 m 2 /s. Continuous contours denote pos¬ 
itive values of the streamfunction, whereas dashed contours 
denote negative ones. The resemblance to Figure 1 is evident. 

FIG. 3. Scatter plot of mean potential vorticity versus 
mean streamfunction obtained from numerical simulation. A 
gross linear tendency appears that resembles the linear ten¬ 
dency observed in inviscid calculations. Each symbol is ob¬ 
tained from a different position in the simulation domain. 
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FIG. 4. Comparison, as a function of the radial wavenum¬ 
ber index Lctk/2n , of the power spectra of (V 2 (ip) + h')/c 
(dash-dot line), (ip) (dashed line), both obtained from nu¬ 
merical simulation, and the power spectra of the solution (ip) 
of Eq. (fjsf) (solid line). 
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